clear all;close all;clc;
N=64;
%=======================================================
% index calcuting
%=======================================================
index = bin2dec(fliplr(dec2bin(0:N-1)));
%=======================================================
% Wnr calcuting
%=======================================================
for r = 0:(N/2-1)
    Wnr_factor  = cos(pi/(N/2)*r) - j*sin(pi/(N/2)*r) ;
    Wnr_integer = floor(Wnr_factor * 2^13) ;
    
    if (real(Wnr_integer)<0) 
        Wnr_real    = real(Wnr_integer) + 2^16 ;
    else
        Wnr_real    = real(Wnr_integer) ;
    end
    if (imag(Wnr_integer)<0) 
        Wnr_imag    = imag(Wnr_integer) + 2^16 ; 
    else
        Wnr_imag    = imag(Wnr_integer);
    end
    Wnr_real_string(r+1,:)  =  dec2hex(Wnr_real,4);
    Wnr_imag_string(r+1,:)  =  dec2hex(Wnr_imag,4);
end
%=======================================================
% fft calcuting
%=======================================================
x_re=1:1:64;
x_im=1:2:128;
fft2 = fft(x_re+x_im*1i)

